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We investigate the propagation of electromagnetic waves in transverse magnetic (TM) mode 
through the structure of materials interface that have permittivity or permeability profile graded 
positive-negative using asymptotic iteration method (AIM). As the optical character of materials, the 
permittivity and the permeability profiles have been designed from constant or hyperbolic functions. 

In this work we show the approximate solution of magnetic field distribution and the eight models 
of wave vector of materials or interface positive-negative gradation. 


I. INTRODUCTION 

A new class of artificial composite material is negative refractive index metamaterials (NRM), which is also called 
left-handed metamaterials (LHM), has attracted the interest of many scientists for more than a decade. Veselagoi who 
first published a review of theoretical NRM for more than four decades ago. Then Pendry translate these concepts 
into practical applications^^— which then attracted great attention on this topic. Experimental confirmation is given 
by Smith et al^ and Shelby et alii. 

The concept of NRM is quite simple. The function of the material is more likely to be determined by the structure 
rather than by chemical composition. This is a new class of materials that goes beyond conventional materials because 
it has the effect of which has never been observed beforeii. 

NRM has structure in the order of subwavelengths and is able to provide negative refractive index in a certain 
wavelength range. The artificial structure contains negative permeability obtained from the double split-ring resonator 
and negative permittivity obtained from nano-wires^. The initial idea of Veselago derived from the dispersion equation 
which expresses the interaction of electromagnetic waves with a dielectric. Dielectric response to the presence of 
an electromagnetic wave is expressed by the basic characteristics of the two quantities, i.e. permittivity, e, and 
permeability, /i. That is because only these two quantities that appear in the dispersion equation^. 

Using a diagram of e — /i, the material properties can be grouped into four classes i.e—: (i) Ordinary materials 
(e > 0, /r > 0); (ii) Materials with negative permittivity and positive permeability (e < 0, /x > 0); (iii) NRM 
(e < 0, /X < 0) and (iv) Materials with positive permittivity and negative permeability (e > 0, /x < 0). 

Analysis of the wave propagation through the NRM structures most commonly done with modeling and numerical 
simulation, in particular using the finite difference time domain method (FDTD). FDTD method has been widely 
accepted as one of the popular numerical methods in computational electromagnetio^. At present, the method used 
to design and to obtain insight into the physical characteristics of electromagnetic NRMi^. NRM structures with 
graded refractive index interface negative-positive has been solved analytically by M. Dalarsson for the special case 
in which its index profile is a linear or exponential functioni^. 

Actually, the real structure of any material contains positive and negative refractive index and at the same time tend 
to have a gradation profile!^. NRM gradation index has been widely studied. Ramakrishna describes an NRM lenses 
composed of media with a gradient index^. Calculation of the transmittance at an exponential gradation structure 
NRM has been given by N. Dalarssoiii^. Analytical solutions obtained from hypergeometry differential equations 
compared with the numerical solution using the transfer matrix method (TMM) seem to have a good pattern match. 
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Electromagnetic wave equation for inhomogeneous medium can be converted into a second order homogeneous 
ordinary differential equation which is linear in one dimensional case. Thus, the analytical solution of Maxwell’s 
equations can be solved in this domain. Homogeneous second order ordinary differential equations appeared widely in 
the literature and there are several ways to obtain an exact solution such as using methods of hypergeometry, super- 
symmetry quantum mechanics^, Nikiforov-Uvarov (NU)^, Romanovski polynomialsi^Ti^ and which has recently been 
used is asymptotic iteration method (AIM) to solve the eigenvalue problem^ and electromagnetic wave propagation 
in inhomogeneous ordinary materiali^. 

Exact solution of held distribution (eigenfunction) and energy (eigenvalue) can be obtained using the AIM if the 
Hamiltonian of the system is shaped like a harmonic oscillator—. In this paper we try to apply the AIM in the case 
of inhomogeneous materials where permittivity or permeability has a positive-negative gradation prohle. Here we 
restrict the calculation just on the approximate solution. 


II. METHOD 

In this section we present mathematical background to manage calculations of light propagation in the interface of 
NRM with positive-negative gradation. We will discuss this study in three subsections: (i) we manipulate Maxwell’s 
equations to construct the wave equation in inhomogeneous materials; (ii) we apply AIM approach to obtain solutions 
of second order homogeneous differential equation; (iii) we apply AIM to obtain approximate solution to the eight 
models of wave vector of materials. 


A. Field Equation 


To simplify the calculations, we set the xz-plane as the plane of ray incident, so d/dy = 0. In addition, the optical 
characteristics of e{x) and /i(a;) are considered as functions of material thickness, x. Using Maxwell’s equations and 
considering the geometry of the material, we obtain a second order differential equation of electromagnetic waves 
in inhomogeneous materials on mode Transverse Magnetic (TM). Actually, in two-dimensional study, there are two 
independent modes i.e. TM mode, the group of fields: {E^, Ey, i?z}, {Ey, Ez, Efj;} , {E^, and Transverse 

Electric (TE) mode, the group of fields: {Elx, Ely, E^}, {Ely, EIz, Ex}, {EI^, Hx, Ey}. 

In principle, these two modes (TM and TE) are similar to each other, so we only need to do one calculation, e.g. 
the TM mode. Now, let’s start by defining the form of plane waves which is a special solution of the wave equation 

H = H{x) (1) 

and from Maxwell’s equations is known that 


iojy,{x) E[ = —V X E] iuje{x) E = x E{ (2) 

Consider the group of fields in the TM mode which we use i.e. [E^, Ex, Hy}. Erom eq.dl]) and eq.© we obtain 


d^Hjx) 

dx^ 


e'{x) dE{{x) 
e{x) dx 


+ y,{x)e{x) — k'^]E[{x) = 0 


( 3 ) 


where in eq. m above, the prime symbol (') indicates the operation of differentiation with respect to x. This second 
order differential equation does not always have the exact solution to a number of functions s{x) and fJ-{x) . In the 
next subsection, we will introduce the approach that is needed to achieve that goal. 


B. Asymptotic Iteration Method 

Asymptotic iteration method aims to obtain the exact solution of a second order differential equation which has 
the following form 

Y"{x) - Xo{x)Y'{x) - so{x)Y{x) = 0 (4) 

where Ao(a:) ^ 0 and so(2:) is the differential equation coefficients and they are so well-defined functions which are 
differentiable. To find the general solution, eq.® needs to be differentiable in {m + 1) times and (m -|- 2) times, where 
m = 1, 2,3,.. is the number of iterations. Then one can obtain 

y(™+l)( 2 .) _ X^_^(^x)Y'(x) - Sjn-l{x)Y{x) = 0 


( 5 ) 
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r(™+2)(3;) _ X^{x)Y'{x) - Sm{x)Y{x) = 0 


( 6 ) 


where 


^m(x} — -|- Atti — I (x) Aq ( 2 ^) “1“ — 


( 7 ) 


Sm(x) = s^_i(a;) + so(a;)Am-i(a;) 


( 8 ) 


which are also called recursive relations of eq. The calculation of the ratio of derivatives of (m + 2) and (m + 1) 
gives 


dx 


In 


y(m+l)(^)' 


y(m+2)(3,) 

y(m+l)(2,) 

Am[A"'(a^) + (5m/Am)F(x)] 
\m-l\X'{,x) + {s„,-i/\m-l)y{x)\ 


(9) 


Furthermore, with the introduction of the asymptotic aspects of this iteration method, i.e. by assuming that the value 
of m is large enough, we obtain 


So eq.® can be written as 


^m{x') SYfi — lix) _ . . 

"x 7 ^ y 7 X “ ot[x) 
Am(x) Xm-l(x) 




Substitute eq. © into eq. (HH and using eq. (HQi), we obtain 

y(m+i)( 3 ;.) = CoXm-i{x) exp|y'[a(x) +Ao(x)] dx 


( 10 ) 


( 11 ) 


( 12 ) 


where Co is the constant of integration. Substitute eq. (ini) into eq.® and solve for Y{x), then we obtain the general 
solution eq.® as 


Y (x) = exp 


— J a{x)dx X ^Ci + Co J exp [2a{x) + Xo{x)]dx'^ dx 


(13) 


where Ci is the constant of integration. 

Eigenvalue of energy is obtained from the combination of termination conditions, eq. (HQD and eqs. O-d 


'^m(x)Sjji—l (x) Xm — l(x)-Smix) — 0, 7Tl — 1, 2, 3 (1^) 

In general, if the problem expressed by eq.® has Xo{x) and so(a7) which are forming a constant function, then the 
general solution of eq. ® can directly follow to eq. (nsi). The most frequent case is that Ao(a;) or so(a;) is not a constant 
function. In such cases, the exact solution eq.® requires a special technique that can be found i n^nd^ie ^ 

Here, we will apply the AIM to the case where the electromagnetic optical character changed gradually. Our 
main attention is the solution domain where the positive-negative interface occurs and how to develop appropriate 
approaches. As an overview of approaches that we use here, let us look at the following example. 


Example: 

Optical characteristics of vacuum space is defined as follows 

e(x) = eo; n{x) = /io (15) 

with £o{x) and are respectively the permittivity and permeability constants of vacuum space, and = 1/c^. 
Substitute ea. (IT5|) into eq.®, we obtain 


d^Hjx) 

dx'^ 


+ (wVc' 


k^)H{x) = 0 


(16) 
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Furthermore, if we let 

H(x) = e~^ (17) 

where ipix) needs to be determined using an iterative procedure. Substitute ea. (ll7l) into ea. ifTHll gives 

+ = o (18) 

ax‘‘ ax 

with 

Ao(a:) = 2a:; so(a^) = — 1 + (19) 

In the solution domain close to zero, eq. m can be approximated by MacLaurin series. The approach is carried by 
taking the first two terms of the MacLaurin series, giving 

Ao(a;) = Ao(0) + Ao(0)a: + ^°^x^ + a:^ + ... 

= 0 + 2a; + 0+ ... 

= 2a: 


so(a:) = so( 0 ) + so( 0 )a:+^^pa:^ + ^2^a;3_p 
= (w^/c^ — k^) — 1 + 0 X a; + —x^ + 0 + ... 

= (u?!(? - k^)-\ (20) 


Substitute ea. (l20)) into eq.([T8|), we obtain 

^ , „ 

dx 

at —(5 < X < S, with 6 = 1. Furthermore, using eq. m, it can be obtained 

Arn^m—1 q[(^ /c k 1) 277.]; 777 = 1, 2, 3, 4, 

According to these conditions, the eigenvalues of energy can be obtained as follows 


OJ 


— k — 1 — 277 —^ kji — 




-(277 + 1); 77 = 0,1,2,3,... 


Using ea.(l23l). eq. m can be written as 

cP'ijj{x) d'ijj{x) 


da: 2 


— 2x- 


dx 


2nip{x) = O', 77 = 0,1,2,3,... 


( 21 ) 

( 22 ) 

(23) 

(24) 


Eg. (1241) known as Hermite differential equation which has a solution in a polynomial of order n. The approximate 
solution of eq. dUl) is 


Hr.{x) = e-^^/\-ire^^ 77 = 0,1,2,3,.. (25) 

at —6 < X < S, with d = 1. Using similar techniques, we develop a method to solve the eight models of interface 
positive-negative gradation as follows. 


1. 1st Model 

Define the character of the optical interface gradation as follows 


e(x) = £o[csch{px) — coth(/9a;)] 
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/i(a;) = — ^Q[csch{px) + coth(/9a:)] 

with p is the gradation parameters. If we let u = px, substitute ea. d^ into eq.® gives 

(PH CSch{u) dH ^ I (uj"^ , 2^ tr _ n 

p du ) 

In the solution domain close to zero, using the MacLaurin series for the hrst three terms, then we obtain 

csch{u) = — 


(26) 


(27) 


(28) 


However, ea. (l28l) is not defined at rt = 0, so we need an assumption. If a is the value of u that close to zero then using 
eg. (1^51) it can be obtained csch(a) = l/cr, and csch{—a) = —1/a. Furthermore, at intervals of —a < u < a, eq.(| 
is considered to have a linear function form which connects the points (—cr, —1/a) and (cr, l/cr), i.e. 


/(m) = 


1 


k^]H = 0 


Using (l29)l and letting z = u/ \/2pcr^, eq. (EH) can be written as 

d'^H dH (2a'^ijp 2cr2 
dz^ ^ dz ~^ \ pP p 

at —cr < u < a, with cr = 1. Furthermore, using the eq. o, it can be obtained 

'2cr^a;^ 2a^ 


\ —2n 


\ _ \ _ 

' pc^ p 

According to these conditions, the eigenvalues of energy can be obtained as follows 

"2cr^w^ 


2a 


pc^ 


Using eq. (PI). eq.dnOl) can be rewritten as 




dH 

—-^5— 2z— -h 2nH = 0; n = 0,1,2,3,... 

dz‘‘ dz 

Eg. (1331) is known as Hermite differential equation. The approximate solution of eg. (l27p is 


(29) 


(30) 


(31) 


(32) 


(33) 


Hrp{z) = (-l)"e^ —(e-^ ); n = 0,1, 2,3,... 


(34) 


2. 2nd Model 

Define the character of the optical interface gradation as follows 

e{x) = —eQ[csch{px) + coth(px)] 


p{x) = p.o[csch{px) — coth(pa;)] 
If we let u = px, substitute ea. (l551) into eq.(I3]) gives 

d^H csch{u) dH 1 /uj^ 2 


du^ p^\P ^ ^ 


(35) 


du'^ p 

Compare to eq.([30|), ea. (l36l) differs only the sign in the middle term of eg. (1311 . So, in the same way we can obtain 


(36) 


d^H dH (2a‘^oP‘ 2a , 9 

.. I I 


dz'^ 


-2z- 


dz \ pP 


Uff = 0 


(37) 
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at —fj < u < a, with a = 1. Furthermore, using eq. da, we obtain 

2(7^ 




n—0 


I + 2n 


pc^ p 

According to these conditions, the eigenvalues of energy can be obtained as follows 

^2cr^a;^ 2(7^ 


p 

Using eq. (12S1), eq.(1371) can be written as 

(PH „ dH 




+ 2z—- 2nH = 0; n = 0,1, 2, 3,.. 


(38) 


(39) 


(40) 


dz^ dz 

Eq. (ITO) is known as Differential Equations of Complementary Error Function. The approximate solution of eq. (l36l) is 


H„(z) = A erf c„{z) + B erf c„(-z); n = 0,1, 2, 3,.. 


with A and B are constant. 


poo poo 

rf Cn(z) = ... erf{z) 

J z J z 


dz 


= 2-’"e"^ 


iFi[l/2(n + l);l/2;z2] 2z iEi[l + (l/2)n; 3/2; z^] 


r[l + (l/2)n] 

is an integral n times the Complementary Error Function, and 


r[l/2(n + l)] 


, . a a(a + l)z^ a(a + l)(a + 2) 

(a)fez'= 


_ 

— 


{h)kk\ 


is the confluent hypergeometric function of first kind, whereas r(z) is the gamma function. 


(41) 


(42) 


(43) 


3. 3rd Model 

Define the character of the optical interface gradation as follows 

£{x) = £o; fj.{x) = fio[csch{px) — coth.{px)] 
If we let z = ^px, substitute ea. dmi into eq.® we obtain 


dz^ p^ \c^ / 

In the solution domain close to zero, it can be approached, tanh(z) = z, in order to obtain 

at —15 < z < i5 with 5 = 1. If we let s = j z, and 

g3 g2 

H (s) = e“8 2 “ ipi^s) 

with p{s)'ip{x) needs to be determined using an iterative procedure, we obtain 

dfii . f .dif , . 

_.Xo{s)--so{sH = 0 


(44) 


(45) 


(46) 


(47) 


( 48 ) 
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with 


Ao(s) = 2s — 


(49) 


and 


so(s) = 1 + — 
P 


4 f 


4a;2 


2/3 


fc2 _ ^2 ^ ^3 _ 1^4 

4 


Using the MacLaurin series for the first two terms, ea. (H^ and ea. dSOl) can be respectively reduced to 

Ao(s) = 2s 


(50) 


(51) 


and 


4 




Substitute ea. (l5T]) and ea. (l5^ into ea. (l4^ . we obtain 


(Pip Q dip 




Ip = 0 


Furthermore, using eq. da, it can be obtained 




+ 2n > ; TO = 1,2,3,4, 


According to these conditions, the eigenvalues of energy can be obtained as follows 


1 + 


4 


p^ \ 4a;2 J 

Using eq. ([55]): eq.([53|) can be rewritten as 


P = —2n —>■&„ = — (—^ ^271 + 1 

2 \ pc J 


Pip dip 

—2s—— \-2nip = 0] 71 = 0,1,2,3,... 

as^ as 

Eg. (1551) is known as Hermite differential equation. The approximate solution of ea. (H5|) is 


i7„(s) = e---(-l)"e« 


ds' 


(e-): 


n = 0,1, 2,3,... 


4- 4b Model 

Dehne the character of the optical interface gradation as follows 

£{x) = eo[csc/i(pa;) — coth(pa:)]; p{x) = po 

If we let z = ^px and substitute eg. (1551) into eq.dSj) we obtain 

d^H csch{z) sech(z) dH A f lP , , , 

^ ^ ^ ^ tanh(z) + P]H = Q 


dz"^ p dz p^ \c^ / 

In the solution domain close to zero, using the MacLaurin series for the hrst three terms, we obtain 

tanh(z) = z 


csch{z) sech{z) = 


1 


(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 














However, eg. (1611) is not defined at 2 : = 0, so it needs an assumption. If a is the value of z that close to zero then 
using eq. (ED it can be obtained csch{a)sech{a) = Xja and csch{—a)sech{—a) = —1/a. Furthermore, at intervals 
of —a < z < a, eq. (ED is considered to have a linear function form which is connecting the points (—cr, —1/cr) and 
{a, l/a), i.e. 


f{u) = 


Using eq. (ED and eq.EDi eq.(ED can be rewritten as 


dz^ a^p dz c2 ^ ^ ® 


(62) 


(63) 


at —a < z < a with cr = 1. If we let s = {lup/z, and 

H{s) = e® (64) 

with '0(s) needs to be determined using an iterative procedure, then substitute ea. dMl) into eg. (1551) we obtain 

^4“ “ V(s)^ - so(s)V'= 0 (65) 

ds 

with 

1 (p^c^ 


Ao(s) = 


cr^p 


and 


so(s) 


= — (p- 


2„2 \ 2/3 


fc2 + 


1 / p^c 


2„2\ 2/3 


3 1 4 

S --S . 


\ 4a;^ J 2cr2p 4a;^ 

Using the MacLaurin series for the first two terms, eq. (ED and eg. (1571) respectively can be reduced to 

and 


cr2p V 4a;2 J 

,, 4/^'VV''b.2 


Substitute eq.(ED and eq.(ED into eg. (155D . we obtain 

.2 2/3 


„2„2 \ 2/3 


‘"v.=o 

da^ a^p \ 4^2 J ds p^ \ J 


1/0 

Furthermore, we let r = {jp'c? j^u) sj v2fT^, eq.(ED can be written as 


Then, using eq. (HI, we obtain 


Yt-ir'Y- —kH' = 0 

dr^ dr p 


8(7^ 


^m^m — 1 — IIyj_Q ^ ^ ^ J , TTl — 1, 2, 3, 4, .. 

According to these conditions, the eigenvalues of energy can be obtained as follows 

= —2n kn = ■^\/pn-, n = 0,1, 2, 3,.. 
p 2cr 

Using eq. (ESD, then eq. ED can be written as 

d^U dib 

—- 2r—+ 2n^ = 0; n = 0,1, 2, 3,... 
dr^ dr 

Eq. (1741) known as Hermite differential equation. The approximate solution of eq. (159^ is 


( 66 ) 


(67) 


( 68 ) 


(69) 


(70) 


(71) 


(72) 


(73) 


(74) 


iJ„(r) = eA77^ ). „ = q, 1, 2, 3,... 


(75) 
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5 . 5th Model 


Define the character of the optical interface gradation as follows 

e(x) = Eq; /r(x) = — iJiQ[csch{px) + coth(pa:)] 
If we let u = px, and substitute ea. (l76)) into eq-ljS]), we obtain 

-^ / '^[csch{u) + coth(M)] + k^ \ H = 0 


p2 ^2 

In the solution domain close to zero, using the MacLaurin series for the first three terms, it gives 

2 u 

csch(u) coth(ti) = —h — 
u 2 


( 76 ) 


(77) 


(78) 


However, ea. (l78l) is not defined at w = 0, so it needs an assumption. If a is the value of u that close to zero then 
using eg. ([781 we obtain csch{a) + coth(tT) = (tr^ + 4)/2tT and csch{—a) + coth(—cr) = —(ct^ + 4)/2(t. Furthermore, 
at intervals of —a < u < a, ea. (l78p is considered to have a linear function form which is connecting the points 
[—cr, — (cr^ + 4)/2cr] and [cr, (cr^ + 4)/2cr], i.e. 


(79) 


(I \ 0-^+4 

Using eq. (Ei), then eq. dZZl) can be written as 

at —cr <u< a with cr = 1. If we let /3 = a;^(cr^ + 4) j 2c^ p^ (j"^, z = and 

H{z) = ^ il){z) (81) 

with ipiz) needs to be determined using an iterative procedure, then substitute ea. (l5T]| into eg. (1501) we obtain 


with 


Ao(z) = 2z — z'^ 


(83) 


and 

= 1 + ^2^ - 

Using the MacLaurin series for the first two terms, ea. (l83ll and ea. (l84l) respectively can be reduced to 

Ao(z) = 2z 


and 


^2 

So(z) — I + p2/3p2 


Substitute ea. (l85|) and ea. (l86|) into eg. (1821) . we obtain 

(Pij} dll) 

dz"^ dz 



e 

pp-i^p^ 


= 0 


(84) 


(85) 


( 86 ) 


(87) 
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Then, using eq. o, we obtain 


— l — n^_Q I 1 T 


+ 2n] ; m = 1, 2,3,4, 


p2/3p2 

According to these conditions, the eigenvalues of energy can be obtained as follows 

1 + fj2/3p2 ^ kn= i0^'^p\/2n + 1; n = 0,1,2,3,.. 

Using ea. dMl) . then eq. (EZl) can be written as 

(Pib dtb 

- 2z-^ + 2n'!/'= 0; n = 0,1,2,3,.. 

Eq. (90) known as Hermite differential equation. The approximate solution of eq. (l77l) is 


r3 2 / 2\ 

Hn{z) = e~““(-l)"e^ — [e~^ j ; n = 0,1,2,3,... 


( 88 ) 


(89) 


(90) 


(91) 


6 . 6th Model 

Define the character of the optical interface gradation as follows 

£{x) = —eo[csc/i(/9a;) + coth(pa;)]; ^{x) = fig 

If we let u = px and substitute ea. dMl) into eq.Q gives 

cPH csch{u) dH 1 f w 
du'^ I 


,-T ^ [csch(u) + coth(M)] + k'^ > H 

du P 


(92) 


(93) 


In the solution domain close to zero, using the MacLaurin series for the first three terms, functions of csch{u) and 
csch{u) + coth{u) can be approximated by the MacLaurin series as ea. d^ and eg. (1751) . so that ea. dMl) can be reduced 
to 


drH u dH 1 


+ 4 


u + k^ 


iJ = 0 


(94) 


dv? a'^p du L V 2cr^ 
at —a < u < a with a = 1. Furthermore, if we let (3 = a;^(cr^ + 4)/2c^p^cr^, z 

H{z) = e"'/V(2) (95) 

with ipiz) needs to be determined using an iterative procedure, then substitute ea. (P5|) into ea. dMl) . we obtain 


d^-ib . , . , 

- sq{z)iI) = 0 

dz dz 


with 


and 


^o{z) = 


2p2/3p 


+ z^ 


So(^:) = 


-rz 


j]2/3p2 2a'^fP-l'^p 4 

Using the MacLaurin series for the first two terms, eq. (EZD and eg. (1551) respectively can be reduced to 

^0(2) = ^2^2/3p 


(96) 


(97) 


(98) 


(99) 















and 
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soiz) = 

Substitute ea. dSni) and ea. dlOOl) into ea. dMl) . we obtain 




f]‘2/3p2 


cPtf) 


+ 


z dip k"^ 


dz'^ Cr2/32/3p /32/3p2 

If we let s = z/eq. (llOip can be written as 

d?ip 


ip = Q 


Then, using eq. (H, it can be obtained 


as^ as p 


/2a^ 


— l — bln—0 ^ ^ I ? ^ 


According to these conditions, the eigenvalues of energy can be obtained as follows 

V^. 


2a^ 2 r. 

- r£ ^ 271 —y Kfi ^ 

P cr 


■; n = 0,l,2,3,.. 


Using ea. (|104|) . then ea. (jl02l) can be written as 

d?ip dip 


+ 2s—- 2nip = 0; n = 0,1, 2, 3,.. 

ds^ ds 


( 100 ) 


( 101 ) 


( 102 ) 


(103) 


(104) 


(105) 


Eg. (11051) is known as the Complementary Error Function Differential Equations. The approximate solution of ea. dMl) 
is thus 


l{fl±£) „3 


Hn{s) = e erf Cnis) + B erf Cni-s)]] n = 0,1, 2,3,... 

with erf Cn(s) refers to eq. dUD. 


(106) 


7. 7th Model 


Define the character of the optical interface gradation as follows 

£{x) = £Q[csch[px) — coth(/oa;)] 


p{x) = pQ[csch{px) — coth(pa;)] 


If we let z = \px, and substitute ea. p07l) into eq.([3]), we obtain 


dfH csch{z) sech{z) dH 4 


ra;2 


dz2 


dz 


tanh^(z) — 


H = 0 


(107) 


(108) 


In the solution domain is close to zero, the functions of tanh(z) and csch{z)sech{z) can be approximated by a 
MacLaurin series in the first three terms, and follow the same way of ea. dSH -ea. dS^ . then ea. dlOSI) can be reduced to 


dH 4: f Lo 


d‘^ H .4. ll±-L ^ t UJ o , O 1 

at —a < z < a with cr = I. If we let s = {fktP' j(?(?Y^^ ea. dI09|) can be written as 


(109) 


( 110 ) 
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with 


Ao(s) = 


2wcr2 


and 


So{s) = —k^ — 
pu} 


( 111 ) 


( 112 ) 


Eg. (11111) and eg. (11121) can be reduced using the MacLaurin series for the first two trems, so that eg. piOl) can be 
written as 


(PH 


c dH 2c 


s-- PH = 0 


ds^ 2ujp ds puj 


Furthermore, if we let r = cjAuiP, ea. (jll3|l can be written as 

d^H ^ dH 8P 


Next, using the eo. (IT4l) . we obtain 


j 2 - = 0 

dr^ dr p 


8P 


1 — ffn—0 ^ k 2 tI j ; 772 — 1, 2, 3, 4 

According to these conditions, the eigenvalues of energy can be obtained as follows 

8p 2 * _ 

- k =—2n ^ kn = — Jpn', n = 0,1, 2, 3,... 

p 2a 

Using eg. (113011 . eg. (111411 can be written as 


j 2 u fju 

-2r—+2nH = 0; 77 = 0,1,2,3,... 

dr 


Eg. (11171) is known as Hermite differential equation. The approximate solution of eo. (11081) is 

i7„(r) = ; n = 0,1,2,3,... 


(113) 


(114) 


(115) 


(116) 


(117) 


(118) 


8. 8th Model 

Define the character of the optical interface gradation as follows 

£{x) = —£Q[csch{px) + coth(pa;)] 

p,{x) = — pe)[csch{px) + coth(pa;)] 
If we let u = px, and substitute eg. (jll9ll into eq.([3]) we obtain 

d^H csch{u) dH 1 (P 


dv? p du ' p^ \ 


[csch(7i) + coth(77)] — k > H = 0 


(119) 


( 120 ) 


In the solution domain is close to zero, the functions of csch(u) and csch{u) + coth(M) are reduced in the same way 
as in eo. (l28l) and eq. dUl), so that eo. (11201) can be written as 


d^H u dH 1 
dP a'^p du p^ 


P / (T^ + 4 
~2a^ 


P -¥ 


£7 = 0 


( 121 ) 
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at —(7 < u < a with a = 1. If we let 7 = a;^(cr^ + 4 )^// 3 ^c^( 2 cr^)^, and z = then ea. (ll21l) can be written as 

( 122 ) 


cPH , , ,dH , , ^ ^ 

_ + Ao(z)—-so(z)i? = 0 


with 


and 


= 2 1/2 


(123) 


So( 2 :) = 


p2^1/2 


— z 


(124) 


Ea. (ll23[) and eg. (11241) can be reduced using the MacLaurin series for the first two terms, so that eg. (11221) can be 
written as 


(PH 


dH 


dz^ dz 

Next if we let s = z/p, then ea. (ll25ll can be written into 


i/ = 0 


Next, using eq. da, it can be obtained 


ds^ ds p 


/2(t2 


— l ^ h ‘2 tIJ 5 ^ — 1 ? 2 , 3 , 4 , 

According to these conditions, the eigenvalues of energy can be obtained as follows 

2(7^ n \/Wi 

-= 2n —>• fc„ = ——; n = 0,1, 2, 3,.. 

p 2a 

Using eo. (|128p . ea. (ll26l) can be written as 


d?' H dH 

—^ + 2s—- 2nH = 0; n = 0,1,2,3,.. 

(JjO lLo 


(125) 


(126) 


(127) 


(128) 


(129) 


Eo. (ll29[) is known as the Complementary Error Function Differential Equations. The approximate solution of eg. (11201) 
is thus 


Hn{s) = A erf c„(s) + B erf c„(-s); n = 0,1, 2, 3,... 
with erf c„(s) refers to eg. (H^ . 


(130) 


III. RESULTS AND DISCUSSION 

It has been shown through the eight models that AIM can be used to find approximate solution of electromagnetic 
wave propagation in the positive-negative gradation of NRM. Some materials that are mathematically can be described 
as a function which has a singular point has been approached by a linear function in the domain such that pretty 
close to the singular point. 


IV. CONCLUTION 

This paper has developed an approach for analyzing the distribution of the magnetic field and electromagnetic 
wave-vector in the positive-negative gradation NRM constructed based on the AIM. As a test of this method, has 
been used as an example of vacuum-space calculation. Differential equation approach results from the eight NRM 
models can be solved by Hermite polynomials and the Complementary Error Function. 
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